13 



Absence of self-organized criticality in a random-neighbor version 

of the OFC stick-slip model 

S. T. R. Pinho, C. P. C. Prado, and 0. Kinouchi 

Institute de Fisica 
Universidade de Sao Paulo 
Caixa Postal 66318 
05315-970, Stio Paulo, SP, Brazil 



> 

Q . e-mail: prado@if.usp.br 



(September 23, 1997) 

CN ■ We report some numerical simulations to investigate the existence of a 

self-organized critical (SOC) state in a random-neighbor version of the OFC 
<^ . model for a range of parameters corresponding to a non-conservative case. 

o 

o 

On 

calculations that agree with our numerical conclusions. 

o 

X\ PACS number(s): , 64.60.L, 05.40, 05.70.L 

S": 

'"O ■ Keywords: SOC, Random Processes, Nonequilibrium Thermodynamics. 



In contrast to a recent work, we do not find any evidence of SOC. We use a 
more realistic distribution of energy among sites to perform some analytical 



The first cellular automaton to display SOC (self-organized criticality), the sandpile 
model |l| , required both a bulk conservation law and a lattice structure with open boundary 
conditions, that were thought to be the main ingredients to allow the building of the spatial 
correlations There are some attempts to construct much simpler and more tractable 

mean-field versions of these systems j§-§]. Olami, Feder and Christensen |7j introduced 
a non-conservative cellular automaton (the earthquake model) that displays self-organized 
critical behavior for a certain range of a parameter a related to the breaking of a conservation 
law. A detailed analysis of this model seems to show that the spatial inhomogeneities 
created by the lattice, the particular boundary conditions, and the global driven mechanism, 
are the essential ingredients for the existence of SOC. In a recent publication, Lise and Jensen 
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T0|1 considered a random-neighbor mean-field version of the OFC model with coordination 
q — 4. According to these calculations, it is remarkable that SOC can still be observed in 
non- conservative cases, for 2/9 = a* < a < a c = 1/4. 

In this work we revisit the random-neighbor version of the OFC model, with an arbitrary 
coordination q (that should be of interest in the context of the investigations of SOC in 



artificial neural networks |TIf , in which case the existence of a regular lattice is not a realistic 
assumption). Unlike the conclusions of Lise and Jensen, we do not find any evidence of SOC, 
except in the conservative case. We show (see figure 1) that the use of lattices that are not 
big enough to allow the occurrence of the largest avalanches can lead to the wrong impression 
that the system does display SOC. Of course, this is an important word of caution about 
all numerical evidences of SOC. Within the mean-field approach, we show that the use of a 
simple and slightly more realistic assumption about the form of the stationary probability 
distribution p(E) of the energies of the sites is enough to completely change the conclusions 
of Lise and Jensen. Using this simplified form of p{E), there is just a very small range of 
parameters where SOC might occur in the non-conservative context. 

The random-neighbor version of the OFC model consists of iV sites with an energy 
Ei < E c , for i = 1,...,N, where E c is a threshold value (the energies of the stable sites 
will be labelled by a superscript — ). During a long time scale, there is a continuous build 
up of all energies until the energy of a certain site i reaches the value E c . Then, site i 
becomes unstable (its energy being labelled by a superscript +), and the system relaxes in 
a very short time scale. The unstable site i transfers an energy aE + (a < 1/q) to q random 
neighbors, which may also become unstable and may generate an avalanche that only stops 
when the energies of all sites are again below E c . The equality a = a c = 1/q corresponds 
to the conservative limit. 



To make contact with the work of Lise and Jensen [1(J, we define a "re-scaled" dissipation 
constant, a q = qa/4. According to the same approximate calculations of Lise and Jensen, the 
critical value a* below which the avalanches have a typical size is given by a* = qj (4q + 2). 
If q = 4, we have a\ = 2/9 = 0.222, as obtained by Lise and Jensen. 



Our own simulations, however, indicate that the density distribution of the avalanche 
sizes p(s) converges to a fixed curve that does not depend on the size of the lattice, even for 
values of a q far beyond a*, given by our mean- field calculations (see figures 1 and 2for p (s) 
and the average avalanche size, (s), repectively). Because of practical limitations, it is not 
feasible to simulate systems closer to the conservative limit, but we strongly believe that the 
collapse of p(s) to a universal curve, and the convergence of (s) to typical value, regardless 
of the size of the lattice, for all values of ot q that we were able to use in the simulations, do 
indicate that the system is not in a critical state. 

The results of numerical simulations for p(E), with q — 4, are shown in figure 3 (for 
simplicity we assume from now on that = a). These results clearly show that p(E) is 
not a simple constant. We then decided to use the same approach of Lise and Jensen fL0| , 
but supposing that the distribution p(E) has the (more realistic) form shown in the detail 
of figure 3, where A p is half the width of each peak, A& is the width of the gaps between 
two peaks, and E* = 3A& + 7A p is the maximum value of E for which p(E) ^ 0. We then 
have 

P + (E + )= E -^ = / P (E)dE. (1) 

J p(E) dE P E c -aE+ 


The lower limit of the integral in the numerator, E c — aE + , can belong to any of the intervals 

that define the four peaks of the distribution, to which we assign the indices i = 1,2,3,4. 

Now we have to consider each one of these possibilities (to simplify the notation, we use the 

superscript i to refer to the value of E c — aE + ). The integrals P+(E + ) have the generic form 

■ , i . (i-l)A b E c aE + 

^ ) = l + L_J_t__£ + (2) 



Therefore, the branching rate a is given by 



a = API 



(i-l)Ab _ a(E+) 
7A p 7A p 7A p 



(3) 



Now, we write the mean-value of the energy of an unstable site, 



where 



E+y = J£L ( 4) 



E c 

E p(E) dE 

E-)= E -^1 . (5) 

p(E) dE 

E c -a{E+) 



Therefore, 



E +V_ Ec [7A p + (»-l)A fc ](l-a) y/yl 

I a(2-a) a(2-a) 2«(2-a)' 1 ' 

where 

y i = 4{S c (l-a)-[7A p + (i-l) A b ]} 2 + 4« (2 - a) [ Xi -U{i-l)\ A p A b , (7) 

with Xi = 24, 26, 32, and 42, for i = 1, 2, 3, and 4, respectively. 

To reach a critical state, we have to impose a > 1. Taking a = 1, and using Eq. (^j), 
instead of Lise and Jensen's approximate result, we have 

7A P (2 + a) + 4 (z - 1) A b - 4£ c (1 - a) + 7a£ c ± 2^ = 0. (8) 

For instance, if we take A p = 0.08 and A& = 0.1, the critical branching condition leads to 
values of a* outside of the physical range (that is, a* > 1/4). Therefore, in this particular 
case, it is physically forbidden to assume that o = 1, so there is no self-organized critical 
state. 

In the conservative limit (for A p — * 0, A& — ► aE c ), the four peaks of p(E) tend to four 
delta functions at ( i — 1) a E c . In this case, it is easy to see that a = 1 leads to the only 
possibility a* = a c = 1/4 (we obtain a* > 1/4 for i — 1,2, and 3). It can also be shown that, 
if we consider the limit A& — > and A p — ► E c / 7, which corresponds to the approximation 
of Lise and Jensen, then a* = 2/9. 

In general, from Eq. (|S|), for all values of i, the regions of parameters associated with 
a < 1/4 are determined by 



Ec - —A p - ^A b < 0. (9) 
24 p 21 v ' 

From this inequality, and the relation E c > 7A P + 3 At, (see figure 3), we see that only in a 
very small region of the parameters A p and A& (see figure 4) there are values of a* in the 
physical range (that is, such that < a* < 1/4). In all of those cases, is very small, and 
the shape of p(E) is very close to the constant form used by Lise and Jensen. 

In conclusion, on the basis of a mean- field argument, supplemented by a more realistic 
approximation for the distribution of energies p{E), we give some analytical indications to 
support our own numerical findings that the random-neighbor version of the OFC stick-slip 
model cannot display SOC outside of the conservative limit. At the time we were writing 
this paper, we came to know about two other works P~2| , [T~5[ ] that lead essentially to the same 
conclusions. 
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Figure Captions 

Fig. 1- a) Distribution of the avalanche sizes (number of topples) for q = 4, a = 0.23, 
and different lattice sizes, L. The number of iterations is n = 2 x 10 6 . For L > 400, the 
curves collapse to the same form, which indicates that there is no self-organized critical 
state; b) Using L = 100, a = 0.23, and n = 4x 10 6 , the statistics is similar for different 
coordinations, q — 4, 6, and 10. However, for large q, and a not too large, the redistribution 
of energy is not sufficient to generate larger avalanches. 

Fig. 2- a) Average avalanche size (s) versus lattice size L, for q — 4, and different values 
of a q (smaller and larger than Lise and Jensen's value, a* = 2/9). For a* = 0.23, for 
example, using bigger lattice sizes (L = 600), we can see the exponential behavior of (s) (in 
contrast to Lise and Jensen's simulations). For a q > 0.24, it is necessary to use much bigger 
lattice sizes to see this type of behavior; b) Average avalanche sizes versus lattice size for 
a q = 0.23, and different coordination numbers, q — 2, 4, and 6. For larger values of q, it is 
sufficient to use a small lattice size to see the exponential behavior of (s) . 

Fig 3-a) Distribution of energy per site p(E) versus energy E, for q = 4, and a q = 0.21, 
0.22, and 0.23. The width of the peaks decreases as a q increases; b) Special form of p(E) 
as used in the calculations, where A p is the half-width of the peaks, A b is the width of the 
gaps, a is the value at the peak, and E c > E* = 7A p + 3A 6 . 

Fig. 4-Space of parameters of p(E) in terms of 7 P = A p /E c and 7& = A b /E c . The 
shaded regions correspond to the intersection between a < 1/4 and (7A P + 3Ab) / E c — 
7% + < I. Depending on the value of E c — aE + , we have (a) E c — aE + e [0, A p ]; (b) 
E c - aE + e [A p + A 6 , 3 A p + A b ] ; (c) E c - aE + e [3A P + 2A b , 5 A p + 2 A b \ ; and (d) E c - aE+ 
e [5A p + 3A 6 ,7A p + 3A 6 ]. 
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